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Abstract 

A quantum-mechanical calculation of conductance in an open quantum 
dot is performed in the Landauer-Biittiker formalism using a tight binding 
Hamiltonian with direct Coulomb interaction. The charge distribution 
in the dot is calculated self-consistently as function of a gate potential, 
for various dot-leads couplings. The interaction is active only inside the 
dot, but not in the leads, its strength being an input parameter. Our 
calculations are complementary to the master-equation approach go 
beyond the "orthodox theory", and account for the size, tunneling, and 
interaction effects in quantum dots. 

1 Introduction 

It is known that the electron-electron interaction (EEI) is important in the trans- 
port properties of small systems such as quantum dots (QD's). The standard 
approach is the so-called " orthodox theory" , which neglects the size quanti- 
zation and reduces the quantum mechanics to the charge quantization inside 
a capacitor In this paper we propose a pure quantum- mechanical method, 
based on the Landauer-Biittiker (LB) formalism and the Hartree approximation 
(HA). We consider a discrete, tight-binding model, which allows the tailoring of 
different shapes of the QD and the presence of a magnetic field in the Peierls' 
fashion. The properties are controlled by several inputs: a) the coupling of the 
QD with external leads, tLu; b) the size and shape of the dot which determines 
the electronic spectrum in the absence of the EEI; c) the interaction strength, 
U. The calculations are performed self-consistently in the HA for dots of rect- 
angular shape with four attached leads, accommodating up to 40 electrons. The 
energy spectrum, transmittance, and dot charging are calculated and discussed. 

2 The formalism 

We model the quantum dot as a 2D rectangular plaquette weakly coupled to 
external leads. In the tight- binding picture the Hamiltonian consists of three 
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terms, describing the isolated QD, for the leads, labeled by a, and H^^ 
for the coupling: 

H = H'^ + Y.Ht + Y.Hk''- (1) 

a a 

Since the role of the leads is only to inject and drain the electrons, we neglect 
the EEI in , and include it only in . Meir and Wigreen and Gartner |^ 
have shown, using different approaches, that in this situation the LB formalism 
remains valid. 

In this approach the Fermi level Ep is fixed by the infinite leads, and the 
external gate Vg is simulated by the site energy in H^: 

^'' = E(^^ + f^E^^)i^)(^i+^^ E (2) 

i j^i " <i,i'> 

< rij > is the mean occupation number of the site j S QD, and U is the param- 
eter describing the EEI. We choose to = ^ (i.e. the energy unit is the hopping 
integral in QD) and denote by <i,i' > the nearest-neighbours summation. In- 
stead of using the whole Hamiltonian (1), it is convenient to describe the open 
QD by an effective Hamiltonian with a non-hermitic term which depends on 
energy E and incorporates the influence of the leads: 

Hj^fj =H° + H°^{E - H^y^H^^ =H° + tHl J2 e"''' > (3) 

a 

where is the hopping energy of leads, r — t^o/tL, k is defined by 2^^ cosfc = 
E, and \a) stands for the dot state where the lead a is attached. 

The matrix elements of the Green's function G'^ (£') = {i \ {E—H^ff+iQ)~^ 
j) are calculated numerically using Eqs.(2-3) and the self-consistency condition 

< Uj >=- [ Im GjjE) dE . (4) 

J-oo 

In the LB formalism the conductance matrix gafj is given by the transmittance 
between the leads a and /3: 

2 2 

5a/3 = jT^p =AjTHlsm''k\G+piEF)W a ^ l3 . (5) 

In the weak-coupling limit r ^ I, the transport problem is reduced to a tun- 
neling problem, and the transmittance shows resonances corresponding to the 
energy spectrum of the isolated QD with the many-body effects included. The 
position, width and heights of the peaks depend on the above mentioned pa- 
rameters. 

3 Results and discussion 

Each transmittance peak corresponds to the addition of an extra electron in 
a Coulomb blockade mechanism (for a review see Ref. 2). The position and 
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width of resonances can be correlated to the Hartree spectrum of the isolated 
QD displayed in Fig.l vs. Ep. The spectrum is independent of Ep as long as 
the later is not sufhciently close to an energy level. When Ep approaches the 
n-th eigenvalue of the system with N electrons, and the addition of the + 1-th 
electron becomes possible, the whole spectrum raises with the additional energy 
Eadd{N, N + 1) = En{N + 1) - En{N). The essential characteristic is that the 
Eadd is not supplied step-like but it is gained linearly with increasing Ep in 
an interval AEp = Eadd (in this region, a slope=1.0 is evident in Fig.l). This 
calculation shows that Eadd{N, A^-|-l) depends on the number of electrons N and 
on the interaction strength U. With increasing U, Eadd increases, the plateaus 
of the energy levels shrink, while those of the number of particles increase. 
A new electron is added exactly in the middle of the interval AEp, i.e. at 
Ep = {En{N) + En{N + l))/2, see the diagonal line of Fig.l. These charge- 
degeneracy points correspond to the resonance peaks in Fig. 2b. The Coulomb 
blockade is accompanied by the so-called (elastic) cotunneling effect meaning 
that the transmittance is non-vanishing also in a range about the degeneracy 
point This effect gives rise to some width of the Coulomb oscillations. 
The comparison of Fig.l and Fig. 2b shows that the width at the bottom of a 
conductance peak equals Eadd: indicating that, for small tj^jj, the cotunneling at 
resonance is of Coulomb origin. The number of electrons in the dot is a rather 
smooth function of the gate potential Vg, i.e. the steps of Fig. 2a are not well 
defined even for a weak coupling (the dotted line), in spite of the fact that the 
transmittance shows sharp peaks. One notices that each of these peaks points to 
an ideal step of Fig. 2b, which corresponds to the limit t ^ 0. For t ^ 0.3 — 0.5, 
the tunneling effects become more important: the intra- valley cotunneling Q 
appears and even the resonance width is now determined by tunneling. With 
still increasing the dot-leads coupling the resonances are spoiled, the charge 
quantization is totally lost and, if a strong enough magnetic field is applied, 
typical quantum Hall effects appear. We mention that here t^o is present in all 
powers of the perturbation theory. 

4 Conclusions 

This approach based on Landauer-Biittiker formalism, tight-binding Hamilto- 
nian and Hartree approximation is able to describe the interplay of the size, 
interaction, and tunnelig in QD. Due to the size effects, the addition energy 
(in other words, the inverse of the dot capacitance) depends on the number of 
electrons, and the Coulomb oscillations of the conductance are distributed irreg- 
ularly. Both intra-valley and at-resonance cotunneling can be put into evidence; 
at resonance, a competition between tunneling and interaction is present. In- 
creasing the tunneling parameter, the charge quantization vanishes faster than 
the resolution of the Coulomb peaks of the conductivity ; the peaks are enforced 
by the charge-degeneracy condition. 
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Figure 2: (a) The number of particles vs. the gate potential for J7 = (the 
noninteracting system) and for U ~ 0.5. With solid lines - the limit t ^ 0, with 
dashed lines - the weakly coupled system, r = 0.1. (b) The conductance gi2 in 
the tunneling regime, for U = 0.5, with solid line, with maxima corresponding to 
the steps in the number of particles. The dash-dotted line shows the conductance 
for [/ = 1. 
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